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Abstract 
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1 Introduction 



One of the more challenging problems in non-integrable soliton dynamics is to find a collective 
coordinate model of the dynamics. In some theories there is a moduli space of exact multi- 
soliton static solutions, and the soliton dynamics at low energy is modelled by geodesic motion 



on the moduli space HIT] , [17| . In others there is a nearby theory with exact static solutions 
and a moduli space, and the soliton dynamics is modelled by geodesic motion on the moduli 
space perturbed by small forces [|lj| [16| . In exactly integrable models, like the sine-Gordon 
model, one can exactly compute what happens in a multi-soliton or multi-soliton-antisoliton 
process. There is no energy loss to radiation in these processes, and solitons do not annihilate 
antisolitons in a collision. However, the problem of interest here involves soliton-antisoliton 
annihilation. There is no moduli space of static solutions available; neither is it close to being 
exactly integrable. 

A method of dealing with solitons, or solitons and antisolitons, which exert forces on each 
other, is constrained minimization ||. One introduces collective coordinates by hand, say the 
positions of the solitons and antisolitons, and then seeks the field configuration of minimal 
potential energy for given positions. The reduced model for the soliton dynamics is partly 
defined by this (minimal) potential energy as a function of the collective coordinates. There 
are various disadvantages of this procedure. Firstly, it is rather ad hoc, because there are 
many possible definitions of the collective coordinates. When solitons are close together, their 
positions are not uniquely defined. Different definitions will lead to different reduced models. 
Secondly, the kinetic energy is not well-defined. It is tempting to use the naive kinetic energy 
expression, based on the soliton masses and the collective coordinate velocities, but this is again 
quite misleading when the solitons are close together. We have learnt from examples like BPS 
monopoles, where there is a moduli space of exact static solutions, that the kinetic energy is 
determined by the geometry of the moduli space, and it could not easily be determined if one 
had a rigid a priori view of what the collective coordinates of the monopoles were. 



The method we use here, advocated in |12| and previously developed in the context of 



instanton-antiinstanton physics in Q, exploits the gradient flow equations associated with the 
field theory of interest. The detailed procedure is explained below. Here let us point out that 
gradient flow represents a dissipative dynamics of fields, and hence of solitons and antisolitons. 
But this gradient flow is only used as an auxiliary tool to construct a non-dissipative reduced 
dynamical system, with its own Lagrangian. The gradient flow curves are chosen to map out 
a finite dimensional manifold of field configurations, and the reduced Lagrangian is simply 
the field theory Lagrangian restricted to this manifold. A natural choice for the gradient flow 
curves is the union of all those that descend from a particular unstable, static solution of the 
field equations. This defines the unstable manifold of the solution. Alternatively, for solitons 
that attract each other, one can define the unstable manifold as the union of all the gradient 
flow curves which descend from configurations of very well separated solitons. 

The unstable manifold acquires a metric and potential from the field configuration space 
of which it is a submanifold. These define a Lagrangian for dynamical motion on the unsta- 
ble manifold, leading to an Euler-Lagrange equation of the form (|3.2|) below. This reduced 
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dynamical system is the collective coordinate model for the soliton dynamics. Any set of co- 
ordinates on the unstable manifold can be regarded as collective coordinates for the solitons, 
although there may be a preferred choice if the solitons are well-separated. A solution of 
the Euler-Lagrange equation on the unstable manifold is not, in general, identical with any 
solution of the original field equations, but it may provide an approximation. 

The reduced system should be particularly useful if the potential energy function does not 
vary much over the unstable manifold, that is, if the static forces between solitons are weak 
(but nonvanishing) , so large kinetic energies are not generated. The unstable manifold is then 
an almost flat valley in the field configuration landscape, and the motion on it at modest 
speeds is expected to be a good approximation to the true field dynamics. This is because 
motion orthogonal to the valley, i.e. radiation, is not readily generated. However, this has not 
really been demonstrated analytically, nor tested numerically. In the example to be considered 
here, the potential energy varies substantially, but we shall show that the reduced dynamics 
is still useful. 

The gradient flow in field theories with solitons has been investigated in detail in some 
cases. Demoulini and Stuart have studied the gradient flow in the two-dimensional abelian 
Higgs model |J, whose solitons are magnetic flux vortices, proving the global existence of 
gradient flow curves, and obtaining some insight into the unstable manifolds. For N type II 
vortices, which repel, the unstable manifold of the static axisymmetric solution representing N 
coincident vortices is particularly interesting. The reduced dynamics on this manifold should 
be a good model for iV-vortex motion, but this has not been investigated. Waindzoch and 
Wambach have studied the gradient flow in the two-Skyrmion sector of Skyrme's soliton model 



of nucleons [18], the geometry of which was discussed in ||, and have numerically (partially) 
constructed an unstable manifold which should model the two-Skyrmion dynamics. Also, 
there has been a detailed study of the gradient flow equation in 4 theory with defined 
on a finite spatial interval [fF|. It has been shown that the unstable manifold of the simplest 
unstable solution on the interval, namely = 0, completed with the unstable manifolds of 
lower energy solutions, is a global attractor for the gradient flow. However, there has been 
no comparison made between the solutions of the 4 theory field equation and the dynamical 
motion on the unstable manifold. This example of 4 theory on an interval shows that an 
unstable manifold can have a non-smooth boundary (see Fig. 8.5 in ref. J7|]). For an unstable 
manifold to be useful for modelling multi-soliton dynamics, and especially quantized soliton 
dynamics, it is important that the manifold has no boundary at all. The vortex and Skyrmion 
examples appear to satisfy this requirement. Mathematically, it is not at all clear in what 
circumstances an infinite-dimensional gradient flow system has a smooth, finite-dimensional 
global attractor without boundary. 

Here we apply the gradient flow method to a particular case of kink-antikink dynamics in 
4 theory defined on an infinite line. The reduced system is one-dimensional and describes a 
process by which the kink and antikink annihilate, exciting a particular mode of the field. In 
the reduced system, this process then reverses and the kink-antikink pair is recreated. Since 
we are interested in how the kink and antikink annihilate into radiation, we perform a further 
analysis of how the reduced system couples to the radiation modes, and estimate the rate at 
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which the energy of the reduced system is transfered to radiation. Fortunately, this rate is 
relatively slow, and we find a good agreement between the complete field dynamics and the 
dynamics of the reduced system coupled weakly to radiation. 



2 </> 4 Theory 

The model we consider is 4 theory on a line. This model is a nonlinear, Lorentz invariant, 
scalar field theory which is not integrable [HJ. The field 4>{x,t) has kinetic energy 



oo 1 . 
- 1 2 



T = / -V dx (2.1) 

-oo / 



and potential energy 



oo 1 1 

- 1 12 i L /i2 i\2\ 



V = J_J-^ + ^-iy)dx. (2.2) 

The Lagrangian is L = T — V, and this leads to the field equation 

- 0" + 20 (0 2 - 1) = 0. (2.3) 

There are two vacua, = ±1. There are also two types of solitons, called kinks and antikinks. 
A kink interpolates between —1 and 1 as x increases, and an antikink does the reverse. Eq.( |2.3| ) 
has static solutions 

4>k{x — a) = tanh(x — a) and 0^(x — b) = — tanh(x — b), (2.4) 

where <pK and <px are t ne kink and the antikink respectively. They have potential energy 
(rest mass) |. The parameter a (b) denotes the centre of the kink (antikink). These solutions 
can be Lorentz boosted to an arbitrary speed v < c (where the speed of light c = 1 in our 
units). There are non-static field configurations with several kinks and antikinks, but kinks 
and antikinks must alternate along the line. A neighbouring kink and antikink can annihilate 
into oscillations of the field - interpreted as radiation. Also, in suitable circumstances, a 
kink-antikink pair can be produced. 

There has been a substantial investigation of kink-antikink dynamics (see, e.g. ref. ||). If 
a kink and antikink approach from infinity, starting with any speed much less than the speed 
of light, they will annihilate into radiation. In higher speed collisions, what happens depends 
sensitively on the speed. 

We have investigated a different process, which appears to be better suited to a collective 
coordinate analysis. This is the symmetric motion of a configuration of two kinks and an 
antikink, which we denote KKK. We consider fields which possess the reflection symmetry 

0( x ) = -0(-x) (2.5) 

at all times, and for which — >• 1 as x — > oo. This allows, among other things, for there to 
be kinks centred at ±ct, and an antikink at the origin. A field of this type can evolve into a 
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single kink plus radiation, but the radiation is not emitted very rapidly, as we shall show. In 
the case of pure KK annihilation, all the energy goes into radiation [] 

An interesting and important feature of 4 kink dynamics is that a single kink has a 
normalizable discrete vibrational mode. The discrete mode deforms a static kink at the origin 
to the form 

sinhf x ) 

cj)(x) = <j> K (x) + A r) D {x) with r) D {x) = — -frL , (2.6) 

cosh [x) 

where, in the linearized approximation, the amplitude A oscillates with frequency cu = y/3. 
In comparison, the continuum radiation modes have frequencies uj > 2. Fields of the form 
( |2.6| ) have the reflection symmetry ( |2.5| ). When A ~ —2, the field ( |2.6|) looks rather like a 
KKK configuration with the kinks close together. In fact it is possible to smoothly interpolate 
between a KKK configuration with the kinks well-separated, and a single kink, via fields of 
the approximate form ( [2.6| ). 



3 Gradient Flow 

We recall that a natural (finite-dimensional) Lagrangian system on a configuration space C 
has a Lagrangian of the form 

L(j, y) = T(y, y) - V(y) , where T(y, y) = ^gMff ■ (3-1) 

Here, y l are arbitrary coordinates on C, and gij{y) can be interpreted as a Riemannian metric 
on C. The Euler-Lagrange equation is 



whereas the gradient flow equation is 



dV 

9ljyJ = ~W ' (3 ' 3) 



The gradient flow is orthogonal to the contours of the potential V - the notion of orthogonality 
requires a metric - and in the direction of decreasing V. If yo is a saddle point of V, then the 
unstable manifold of yo is defined as the union of all the gradient flow curves that descend 
from y . If the configuration space is Euclidean, and y l are Cartesian coordinates, then the 
gradient flow equation is obtained from the Euler-Lagrange equation by replacing second time 
derivatives by first time derivatives. 

(/> 4 field theory can be thought of as an infinite-dimensional Lagrangian system whose 
configuration space C is the space of fields {(fi(x)} at a given time. The simple form of the 



1 This is also a rather slow process because an intermediate, long-lived breather state is produced || . 
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kinetic energy expression (|2.1| ) implies that C is Euclidean, with the Riemannian distance 
between fields 0(x) and 0(x) + 5<fi(x) being the square root of 

(5<f)(x)) 2 dx . (3.4) 



The gradient flow equation is therefore obtained from eq.( |2.3| ) simply by replacing by 0, 
giving 

= 0" -2(f) (0 2 -l). (3.5) 

Here, we solve this gradient flow equation to construct a collective coordinate model for the 
KKK system, with the imposed reflection symmetry (|2.5| ). Because of this symmetry, the 
KKK system is modelled by a one-dimensional unstable manifold. It consists of the field 
configurations 0(x), parametrized by — oo <t< oo, occurring in the solution of the gradient 
flow equation ( |3.5| ) which descends from a KKK configuration, with kinks at ±oo and an 
antikink at the origin, to a single kink at the origin. This solution approaches the single 
kink tangent to its discrete mode. The collective coordinate can be identified either with 
the separation of the kinks from the antikink, while this is large, or with the discrete mode 
amplitude, while this has a modest or small negative value. 

The unstable manifold is infinitely long. At the top end, where there are kinks at infinite 
separation, the potential energy is 4 (three times the kink rest mass). At the bottom, where 
there is a single kink, the potential energy is |. Most of the unstable manifold consists of well- 
separated KKK configurations, with kinks at ±a and an antikink at the origin. We cannot 
write down the exact form of the fields satisfying ( |3.5| ), but for a » 1 an exponentially 
accurate approximation is given by the product form 

= -4>k{x + a) 4>k(x) 4>k{x - a) , (3.6) 



with a a function of time. To a sufficiently good approximation, the Lagrangian for such fields 



is fa 2 + 32e~ 2a - 4, since there are two kinks moving with velocities ±d, and the interaction 



potential of a KK pair separated by a is — 16e 2a if a is large JTOf- The gradient flow equation 
therefore reduces to 

a = -24e" 2a (3.7) 
whose integral is e 2a = — 48t + const. We have verified this by taking the initial data (|3.6|) 



with a = 5, and evolving it numerically with the gradient flow equation ( |3.5|) until the three 
zeros of are at x = ±4 and x = 0. The time this takes is very close to ^(e 10 — e 8 ) ~ 397, 
which is a very long time compared with the period of the discrete mode oscillation. 



For a < 4, the expression ( pT6| ) fails to be a good approximation for the fields on the 
unstable manifold, so we have integrated ( |3.5D numerically, starting at t — with initial data 
( |3.6| ) with a = 4. We find that after a finite time the KKK configuration annihilates into 
a deformed single kink. An exact annihilation time can be defined as the time when the 
three zeros of coalesce into one. Equivalently, this is when 0' is zero at the origin. (0' is 
negative before this time and positive afterwards.) The annihilation time, for our initial data, 
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is t — 59.64. By integrating (|3.7Q and finding when a = 0, one may estimate the annihilation 
time to be t = ^(e 8 - 1) ~ 62. 

As t increases further, the field configuration rapidly approaches the single kink, and it does 
so tangent to the discrete mode, since this is the lowest frequency mode of vibration around 
the kink. So the field becomes well approximated by the form ( |2.6| ), with an amplitude A(t) 
tending to zero. Even at the annihilation time, the field (|2.6|) with A = —1 is a reasonable 
approximation, since for this value of A, <p' vanishes at the origin. The asymptotic behaviour 
of A(t) is given by calculating from ( |2.1| ) and ( |2.2j ) the kinetic and potential energies for fields 
of the form Q2.61 ). These are, respectively, 

T(A) = l -A 2 (3.8) 

and 



giving a gradient flow equation 



V(A) = ^ + A 2 + ^A 3 + ^A\ (3.9) 

I A = _2A - —A 2 - —A 3 . (3.10) 
3 8 35 V ; 

As t — > oo, A decays like e~ 3t , the coefficient being the square of the discrete mode frequency. 

Our numerical solution of the gradient flow equation was obtained using a predictor- 
corrector finite difference scheme. Since this is an implicit scheme, there is no stability restric- 
tion on At/ Ax 2 (see |1[). The resulting tridiagonal linear system is solved explicitly using the 
Thomas algorithm. Our space-time domain is defined on —10 < x < 10 and < t < 65 in 
space and time steps of 0.02 and 0.005 respectively. Fig. [I] shows the field <p a t the initial 
time and various subsequent times. The annihilation process occurs at t ~ 60. Fig. |2] shows 
the potential energy as a function of time. The rapid decrease occurs close to the annihilation 
time. 

The unstable manifold stops abruptly at the single kink. It is desirable to define a valley 
which smoothly continues it. The valley continuation needs to approach the single kink tangent 
to the discrete mode from the opposite direction. The best continuation we can think of is the 
set of fields ( [2.6|) with A > 0. With this ansatz, the potential energy in the valley, and its first, 
second and third derivatives, are all continuous at the single kink, but the fourth derivative 
is probably discontinuous. The positive value of A for which the potential energy expression 
(fO) has value 4 is 1.3. 



Any field <p(x), with the reflection symmetry and boundary conditions we are assuming, 
has a decomposition 

<p(x) = <p K (x) + A 7] D (x) + 7](x) (3.11) 

where rj(x) is a superposition of continuum modes orthogonal to the discrete mode t]d{x). We 
can calculate the amplitude of the discrete mode A by projection 

3 r°° 

A = 77 / (<p(x) - <j> K {x)) VD (x)dx . (3.12) 
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In the valley, A varies between —3ir/2 and oo. The first number is calculated using ([3.12j) with 
4>{x) a KKK configuration with kinks at infinite separation, and its finite value implies that 
the tangent to the unstable manifold is asymptotically orthogonal to the discrete mode (both 
metrically, and in the usual sense of integration). The KKK half of the valley is therefore 
curved, turning a right angle along its length. The half of the valley defined using the ansatz 
), with A positive, is straight. 



Along the valley, the arc length s is the only intrinsic geometrical quantity. Suppose s = 
at the bottom, with s positive in the KKK half of the valley. Let V(s) be the potential energy 
in the valley, as a function of arc length. The intrinsic gradient flow equation is 

i = -£ (3.13) 

and in the KKK half of the valley this must give the same dependence of V with time, as 
the field gradient flow equation. Eq. ( |3.13| ) implies 



d^ 
~dt 




(3.14) 



We have taken our numerically obtained V(t) and integrated y — V to find s(t), and hence 
determined s and V(s) along the curved half of the valley. The kinetic energy expression 
( |3.8| ) for a field of the form ( |2.6| ) implies that for small s, we may identify s = —J~^A. 



For large positive s we may identify ds = y|<ia. Numerically we estimate that for large a, 

s = yfa — 1.18. The arc length along the straight half of the valley is exactly s = —^J~^A, 
and using this relation we can easily convert the potential function ( |3.9| ) into a function of s. 
Fig. |H shows the potential function V(s) in both halves of the valley. If we use the discrete 
mode ansatz ( |2.6| ) also for s > we get the dashed curve in the Figure, which is not a bad 
approximation for \A\ < 2. Note the asymmetry of the potential, due to the A 3 term in V. 



4 Dynamics 

The one- dimensional dynamical system 

L = \s 2 - V(s) , (4.1) 

with V(s) as shown in Fig. || is our collective coordinate model for the dynamics of the KKK 
system. It predicts that for energies | < E < 4, there is a nonlinear oscillatory behaviour. For 
E > 4 it predicts that two kinks can approach from infinity and annihilate an antikink at the 
origin, with the kinetic energy being captured by the discrete mode of the kink that remains; 
then the process reverses and KKK reform. The model also predicts that if the initial data 
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is of the form (|2.6|) with A = and \A\ > 2v^2, then a KKK configuration will form, and the 
kinks move off to infinity. (The more direct process is with A initially negative.) 

We do not expect the field dynamics to be exactly along the part of the valley defined by 
the gradient flow, since, as we have shown, the valley is extrinsically curved. However, we 
expect the motion along the valley to be a guide, rather like the bottom of a bobsleigh run 
is a guide to the bobsleigh trajectories. As the valley turns, so the dynamical motion must 
climb up the side of the valley to be forced round the corner. We also expect some vibrational 
motion orthogonal to the valley to be generated. 

We have numerically studied the field dynamics of the KKK system, given by eq.( |2.3| ). 
The integration is performed using a second-order three level implicit formula. As for the 
gradient flow case, this scheme is stable for any value of At/ Ax provided that the three level 
parameter is greater than 1/4 (see ]TJ). Our space-time domain is defined on —20 < x < 20 
and < t < 60 in space and time steps of 0.02 and 0.015 respectively. We have considered 
two kinds of initial data. 



First, we have taken the initial data (|3.6| ) with a = 4 and with = 0. The initial potential 
energy is very close to 4. We observe (see Fig. |]) that the two kinks approach the antikink 
at the origin and annihilate, producing a kink with an excited discrete mode. The system 
oscillates back to KKK, but the separation is reduced. The motion continues, essentially as 
a large amplitude nonlinear oscillation of the discrete mode. Energy is slowly transferred to 
radiation modes, which is emitted symmetrically and escapes to infinity. 

Fig. |5] shows the potential energy and total energy as a function of time. Total energy is 
conserved until t ~ 40, which is when radiation first arrives at the boundary of the simulation 
interval. Our boundary conditions absorb this radiation. Notice that at t = 18.75, after one 
oscillation of the field, the potential energy is about 1.38 which is very close to |, the energy of 
the static kink, so the field dynamics passes very close to the bottom of the valley at that time. 
Notice also the very long time during which the system oscillates with frequency just less than 
V^, and the slow production of radiation. In one period of oscillation the potential energy 
has two minima but these have unequal spacing in time. We have estimated theoretically the 
rate of energy loss to radiation from the oscillating discrete mode, and this agrees well with 
what we see numerically (see Appendix A). 



Fig. |B| shows A(t), the amplitude of the discrete mode calculated using eq. fl3.12p , for our 



dynamical field (j)(x,t). If the KKK system moved exactly in the valley, with total energy 
less than 4, A would oscillate without loss of amplitude. We see in Fig. |] that A oscillates 
with slowly decreasing amplitude. The asymmetry of the oscillation can be understood from 
the asymmetry of the potential about the bottom of the valley. Note that the first time A 
is positive, it reaches 1.15. This is close to the value 1.3 predicted for motion in the valley 
with total energy 4. Fig. |7| shows the field (f) obtained dynamically, at the moment when A as 
defined by ( |3.12 ) is 1.15, compared to the ansatz ( |2.6| ) with A = 1.15. 



Second, we have considered the initial data = 4>k(x), (ft = —Cr] D (x), that is, a single kink 
with its discrete mode excited. This motion starts at the bottom of the valley and tangent 
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to the valley bottom, and although the valley is curved, we expect the motion to be guided 
round the curve, so that a KKK configuration will be produced if C is sufficiently large. 

If C — 2y/2, then the initial kinetic energy is |, which is just sufficient to produce a kink- 
antikink pair. In practice, some energy goes into radiation, and this value of C is not large 
enough. The critical value for kink-antikink pair production is C ~ 3.06. Here, two kinks 
reach the boundary of our region —20 < x < 20, leaving an antikink at the origin, and since 
the forces between the kinks and the antikink are negligible at these separations, the kinks 
would presumably travel to infinity. Since 3.06 is not much greater than 2\/2 , we see that 
exciting the discrete mode of a single kink is an efficient mechanism for producing a KK pair. 
As C increases beyond 3.06, so the outgoing kinks have higher speed, and there is also more 
radiation. Energetically, two KK pairs could be produced when C = 4 but in practice, more 
energy is needed, and the production first occurs when C = 4.71. 

Figs. | and | show the field <j>(x, t) for -20 < x < 20 and < t < 60 with C = 3.06. One 
can see how the discrete mode evolves into the KKK system. The small amount of associated 
radiation is visible too. Fig. [lO] shows the creation of two KK pairs when C = 4.71. 



5 Conclusion 

The KKK system with reflection symmetry is well-described by a one- dimensional reduced 
dynamical system with an asymmetric potential. This is quite surprising because annihilation 
of a kink and antikink can and does occur, releasing a large amount of energy. However, 
most of this energy is converted to the oscillation of the discrete mode of the remaining kink. 
The collective coordinate of the one-dimensional system can therefore be interpreted as the 
separation of the kinks from the antikink, or as the amplitude of the discrete mode, in different 
regions. We have defined the reduced system on one side of the minimum of the potential by 
gradient flow, on the other by an ansatz just involving the discrete mode, but perhaps a better 
ansatz is possible. The reduced dynamics is a good approximation to the full field dynamics, 
because the discrete mode couples quite weakly to radiation modes. 

We conjecture that in other field theories, multi-soliton and soliton-antisoliton dynamics 
can be similarly reduced to finite-dimensional systems using gradient flow, and that discrete 
vibrational modes of the solitons may again play an important role. Soliton-antisoliton pair 
production may be associated with large amplitude excitations of the discrete modes in these 
other field theories. Perhaps soliton-antisoliton pair production at high temperatures can be 
catalysed by the presence of solitons, if these solitons have suitable discrete vibrational modes. 
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Appendix A 



Radiation from an oscillating kink 

In the linear approximation, the discrete mode of the kink oscillates harmonically with fre- 
quency v3- However, when the nonlinearity of the fields is taken into account, this mode 
couples to the continuum radiation modes and eventually all its energy is radiated away. We 
can estimate the rate at which this happens as follows. 

Let us write the field as 

<j)(x, t) = <j) K (x) + A(t) T} D {x) + ri(x, t) (A.l) 

where 4>k{%) — tanh(x) and t]d(x) = sinh(x)/ cosh 2 (x). The function r}(x,t) represents the 
continuum part of the deformation of the kink, and hence is orthogonal to the discrete mode 

r](x,t)r] D (x)dx = 0. (A. 2) 



It is consistent to assume r)(x,t) is odd in x. If we substitute ( |A.1|) into the field equation 
Q2J we find 

(A + 3A)r ]D + 7i-r ] " + (64> 2 K -2)r ] = -6(77 + <j) K )rf D A 2 

-6(77 + 2(p K )r]r] D A 

-2r, 3 D A 3 - 6<W - 2 V 3 (A.3) 

At O(A), the discrete mode oscillates with frequency a/3, and there is no source for 77, so it is 
consistent to set 77 = 0. At 0(A 2 ), there is a source for 77, the term —Q<PkVdA 2 ■ If A is small, 
we may suppose that 77 is 0(A 2 ), and we may neglect terms in ( |A.3|) involving A 3 , rj 2 , rjA, etc. 
The reduced system is 

(A + 3A) VD + 77 - 77" + (6<^ - 2)r/ = -6<f> K r) 2 D A 2 (A.4) 

The source term — QtpKrfbA 2 splits into two parts, a projection onto the discrete mode, and a 
projection orthogonal to this. The projection onto the discrete mode implies an anharmonic 
oscillation of A, but this is not a large effect and the frequency of the oscillation is unchanged to 
lowest order. The projection orthogonal to the discrete mode is the source for 77. Fortunately 
the eigenf unctions of the Schrodinger operator 



_d^_ 

dx 2 



+ (64>k - 2) (A.5) 



are known exactly Jl3| , so the Green's function can be computed and the response of 77 to the 



source can be calculated by explicit integration. 



The truncation ( |A.4j ) is not energy conserving. The discrete mode oscillation creates 



radiation but there is no backreaction on the amplitude of the discrete mode. We may calculate 
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the backreaction by finding the energy carried away by the radiation, and may infer the rate 
of decrease of the amplitude of the discrete mode, using energy conservation. Our calculation 
depends on the rate of energy loss being small during the period 

Let the projection of 4>kVd onto the discrete mode be ar] D , so the orthogonal projection 
is (PkVd ~ ar lD- The constant a is determined by 



4>K(x)r]D(x) 3 dx — a / r]£)(x) 2 dx = 



(A.6) 



implying a — ||. The oscillation of the discrete mode, to 0(A 2 ), therefore obeys 



(A.7) 

and (|3.9| ). The equation 



A + 3A + 6aA 2 = 

which agrees with what would be obtained, at this order, from 
for 1] is 

V~V" + (60x - 2)77 = 6 (arj D - (j) K r} 2 D ) A 2 (A.8) 

Let us now assume that A = A cos( v / 3t), so A 2 = ±A§(cos(2v^t) + 1)- The response of 77 
to the time-independent source is itself time-independent and carries away no energy. The 
important part of 77 is at the single frequency 2 a/3. Let us therefore consider the equation 



r\ — 77" + (6(p 2 K — 2)77 = f(x) exp(iu>t) 
where to > 2. Setting r)(x,t) = rj(x)exp(iut), we obtain 

-77" + (6<&-2 -u 2 )r l = f(x) . 
The homogenous equation ( [A.10Q (with / = 0) has exact solutions |13 



(A.9) 



(A.10) 



Tj q {X) 



(3(j) 2 K (x) — 1 — q 2 — 3iq4>K(x)) exp(+iqx) 
(3(f) 2 K (x) — 1 — q 2 + 3iqcj) K (x)) exp(-iqx) 



(A.ll) 



where q = \Joj 2 — 4. From the form of these solutions as x — > 00, we compute the Wronskian 



W = r] q (x) rf Jx) - r/Jx) rj^ q (x) = -2iq(q 2 + l)(q 2 + 4) . 



(A.12) 



The solution of ( A.10 ) that we want is the one with outgoing radiation. The relevant Green's 
function is 

Y (A.13) 



So the desired solution of (|A.10|) is 

V (x) = -^v~ q (x) f fitKim - ^v q (x) J™ m)v- q (m ■ (a.u) 
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We may use the asymptotic form of (|A.ll ) to obtain the form of the solution ( A.14 ) for large 



x. 



2iq(2 — q 2 — 3iq) 

Thus, for the source —3Al(p K r]1 ) cos(2^/3t), i](x,t) is asymptotically 

/ —3A 2 r°° \ 

V(x, t) = Re L (2 _ g2 °_ 3zg) exp i(cut - qx) j ^ 0^(0^(0^(0^) (A.16) 

where to = 2\/3 and q = 2\p2. It is sufficient here to write the source as proportional to 4>kVdi 
since the integral of f]DVq vanishes. 

Eq.( [A.16|) represents outward moving radiation for large x. There is similar outward radi- 



ation for large negative x. To calculate the integral in (|A.16|) , we use the result 



/oo 
M0vU0(3<Pk(0 -l-q 2 - 3*g<M0) exp(^0d£ 
-oo 



' — oo 

J2 



48 sinh(7rg/2) 
Thus ( |A.16|) reads now 



111(1 (g 2 + 4)(^-2). (A.17) 



V(x, t) = * qiq u % \ ^TT A ° co < ut -W- 6 ) (A.18) 
32 sinh(7rg/2j y 9 + 1 

where 5 is a phase depending on q. We are only interested in the amplitude of the radiation, 
not its phase. For q = 2a/2 the amplitude is 

R= \ n A 2 Q = 0.0453 Al (A. 19) 

sinh(7rv2) 

For a wave of the form r] = Rcos(ut — qx — 5) the average energy flux to the right is ^R 2 uq, 
so the total energy flux away from the oscillating kink is R 2 uq. The rate of energy loss by the 
discrete mode (averaged over a period) is therefore 

— = -(0.0453) 2 4 v / 6 Aq = -0.020 A^ . (A.20) 

Lit' 

On the other hand, the energy of the discrete mode is E = Aq, so the rate of decay of the 
amplitude is given by 



dt 

which can be integrated to give 



a 42 

- -0.020 Aq (A.21) 



1 1 0.020 (t - 1 ) . (A.22) 



A 2 (t) A 2 (0) 
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We can calculate the decay of the discrete mode amplitude due to radiation in another 
way. Consider again the exact equation ( |A.3| ). 

Using the orthogonality condition ( A.2|) , we project the equation ( A.3|) onto the discrete 
mode 

A + u 2 qA = a + a x A + a 2 A 2 + a 3 A 3 (A.23) 

where we have defined 



a 



/OO rOO 

(30x + r])r] 2 r) D dx; ai = -9 / (2cj) K + rf)rrq 2 D d.x 
-oo J — OO 

OO 



a 2 = — 9 / r]7] D dx — 6a; 



—oo 

00 4 - 12 



a 3 = -3 / i] D dx = - 



(A.24) 



35 



Following ref.||, and as in we seek a perturbative solution for (j)(x,t) of the form (|A.1| ) 
with 



A = A^+A^+A^ + --- 
rj = rf x) + rj {2) + r/ 3) H 



and with A having period 27r/u, where 



CO 



(A.25) 



(A.26) 



At the first order, the discrete mode oscillates with the frequency ujq = and there is no 
source term for the continuum modes, so A^ = Aq cos(cuot) and = . 

At the second order, eq. (|A.23|) reduces to 



i (2) + u 2 A ( V = -6aA^ 2 + 2u u {1) A cos{u t) . 



(A.27) 



The requirement that the resonance term (at frequency oj ) be absent implies a/ 1 ) = 0, so we 
obtain exactly eq.( |A.7| ). The latter displays the first anharmonic correction to the discrete 
mode oscillation. Writing A^ = ^Al(cos(2oj t) + 1), the solution of eq.( |A.27|) reads 



A i2) = —A 2 (cos(2u t) - 3) . 



(A.28) 



The correction ( |A.28|) tells us that the main source of the asymmetric oscillation of A in Fig. 
D is the cubic term of the potential (|3.9| ). 

By expanding eq.( [A.3|) up to the second order in r] and A, we obtain 



which is exactly eq. (|A.8|) , therefore t/ 2 -* is given asymptotically by (|A.18|) . 



(A.29) 
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^From eg. ( |A.27| ) we see no r) contribution to the correction, consequently there is no 
energy conservation at this order. To find the modification to the discrete mode oscillation 
caused by radiation, we must go to the next order. The third order part of (|A.23| ) is given by 



A^+oo 2 A^ = (6^i? 2 -6^ + ^)^cos(3^) 

u 4 



a 2 3 

+ (2oo J 2) + QiqR 2 A 2 + 30— A 2 , + -a 3 A 2 )A cos(w t) (A.30) 

(JJ Q I 



where q = 2\/2 and R = R/A 2 , with R given by ( |A.19| ). Note the appearance of complex 
numbers resulting from the coupling of the discrete mode with the imaginary part of the 
continuum (the real part doesn't contribute at this order). The condition for the resonance 
term to be absent is 

= -(SO^ + ^3 + 6^ 2 )^ 
uJq 4 2u 

= -(0.176 + i0M0)A 2 (A.31) 

The second order real correction to the frequency is negative, so the period of the discrete 
mode oscillation is greater than 2tt/ujq. This effect is confirmed by Fig. |6| where we see a 
longer period at earlier times, when the discrete mode amplitude is large. The imaginary 
correction to the frequency implies a decay of the amplitude, the effect we are seeking, which 
in turn means that the real frequency increases with time. The decay rate of the amplitude is 
given by 

^ = -0.010 Al (A.32) 

which is equivalent to (|A.21|) . Fig. [ll] shows the evolution of <f){x,t) calculated numerically 
with initial data (j)(x, t) = 4>k{.%) + T]d{x), so A (0) = 1. The dashed line shows the amplitude 
calculated up to the second-order, and including the frequency correction ( |A.31|) . The decay 



of the amplitude of the discrete mode is thus seen to be well approximated by (A. 22) 
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Figure captions 



Figure 1: The gradient flow annihilation process for the KKK system at time (solid), 59 
(dotted), 59.64 (dashed) and 61 (dashed dotted). 



Figure 2: Time dependence of the potential energy along the gradient flow curve. 



Figure 3: Arc length dependence of the potential energy in the valley. Gradient flow (solid) 
and discrete mode results for s < (the discrete mode result for s > is also shown). 



Figure 4: The full dynamics annihilation process. 



Figure 5: Time dependence of the potential energy (solid) and the total energy (dashed). 



Figure 6: Time dependence of the discrete mode amplitude. 



Figure 7: The field (solid) compared to the discrete mode ansatz at A = 1.15 (dashed). 



Figure 8: Creation process of three kinks for C = 3.06. 



Figure 9: Snapshots of the three kinks creation process at time 10 (solid), 30 (dotted), 60 
(dashed) and 90 (dashed dotted). 



Figure 10: Snapshots of the five kinks creation process at time 10 (solid), 20 (dotted), 40 
(dashed) and 60 (dashed dotted). 



Figure 11: Time dependence of the discrete mode amplitude with Aq(0) = 1 (solid) compared 
to the second order approximation (dashed). 
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Figure 1: The gradient flow annihilation process for the KKK system at time (solid), 59 
(dotted), 59.64 (dashed) and 61 (dashed dotted). 
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Figure 2: Time dependence of the potential energy along the gradient flow curve. 
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Figure 3: Arc length dependence of the potential energy in the valley. Gradient flow (solid) 
and discrete mode results for s < (the discrete mode result for s > is also shown). 




Figure 4: The full dynamics annihilation process. 
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Figure 5: Time dependence of the potential energy (solid) and the total energy (dashed). 
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Figure 9: Snapshots of the three kinks creation process at time 10 (solid), 30 (dotted), 60 
(dashed) and 90 (dashed dotted). 
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Figure 10: Snapshots of the five kinks creation process at time 10 (solid), 20 (dotted), 40 
(dashed) and 60 (dashed dotted). 
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Figure 11: Time dependence of the discrete mode amplitude with A (0) = 1 (solid) compared 
to the second order approximation (dashed). 
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